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Abstract 
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broken BRST symmetry directly translates into a set of Slavnov- Taylor identities. We show that 
these identities guarantee the multiplicative renormalizability of both Gribov-Zwanziger and Refined 
Gribov-Zwanziger theories to all orders. The known property that only two renormalization factors are 
needed is recovered. The non-renormalization theorem of the gluon-ghost-antighost vertex as well as 
the renormalization factor of the Gribov parameter are derived within the linearly broken formulation. 



*capri@ufrrj.br 
t ajgomez@uerj.br 
*msguimaraes@uerj .br 
Vitor@dft.if.uerj.br 
^sorella@uerj .br 

' Work supported by FAPERJ, Fundacjao de Amparo a Pesquisa do Estado do Rio de Janeiro, under the program Cientista 
do Nosso Estado, E-26/101. 578/2010. 
* * dgtedesco@uerj.br 



1 



1 Introduction 



The issue of the BRST symmetry versus the color confinement is one of the highhght of the current 
debate on the infrared behavior of Yang-Mills theories. 

Needless to say, the BRST symmetry is a fundamental tool in order to: i) prove the renormalizabil- 
ity of the theory, ii) identify the physical subspace whose states have positive norm, in) guarantee that 
the restriction to the physical subspace of the scattering operator S is unitary. All this applies to the 
perturbative regime and to nonconfining gauge theories for which the asymptotic fields and their corre- 
sponding elementary particles can be safely introduced. 

Things become much more complicated when the theory under investigation is a confining theory, such as 
pure, i.e. quarkless, Yang-Mills theories. Gluons are not part of the spectrum, due to the nonperturba- 
tive phenomenon of color confinement. The contact with the physical spectrum of the theory is encoded 
in the correlation functions of suitable composite colorless operators built out from the elementary gluon 
field. At present, it is safe to state that we lack a clear understanding of the role played by the BRST 
symmetry in a confining theory. 

In this work we aim at pursuing the investigation on the BRST symmetry within the so called Gribov- 
Zwanziger framework O O |4], which enables us to take into account the nonperturbative effect of the 
Gribov copies by restricting the domain of integration in the functional integral to the first Gribov hori- 
zon, see [5] for a review. Remarkably, this restriction can be implemented within an Euclidean field 
theory framework. The corresponding action, known as the Gribov-Zwanziger action, enjoys the prop- 
erty of being local and renormalizable [3l|4l[6l[7l[8l|9]. As a consequence of the restriction of the domain 
of integration to the Gribov horizon, the gluon two-point correlation function gets deeply modified in the 
infrared region, namely 

where 7 is the Gribov parameter [21 [3l HI [5]. Expression ([1]) displays complex pole, so that it cannot 
describe a physical particle, being suitable for a confining phase. 

The issue of the BRST symmetry within the Gribov-Zwanziger framework has been much studied in 
recent years [TOl [HI HI HSl [HI [15l |16l [TTJ [1], as summarized below. The first observation which can 
be made is that the Gribov-Zwanziger action does not enjoy exact BRST invariance. The usual BRST 
symmetry of the Faddeev-Popov action turns out to be broken by terms proportional to the Gribov pa- 
rameter 7^. Moreover, as the resulting breaking term is of dimension two in the fields, it is a soft breaking 
im I12j. As such, it can be kept under control at the quantum level. This amounts to treat the breaking 
as a composite operator which is introduced into the theory through a set of external fields. This en- 
ables us to write down generalized Slavnov- Taylor identities which guarantee the renormalizability of the 
Gribov-Zwanziger theory to all orders. These generalized Slavnov- Taylor applies as well to the renormal- 
ization of gauge invariant composite operators [18]. In \13\ 115). the softly broken BRST symmetry of the 
Gribov-Zwanziger action has been cast into an exact invariance which is, however, non-local. Further, in 
|16] this non-local invariance has been reformulated as an exact local symmetry by introducing a suitable 
set of additional localizing fields. Though, the resulting BRST operator is not nilpotent. Attempts at in- 
terpreting the breaking of the BRST symmetry as a spontaneous symmetry breaking can be found in [T7] . 

Recently, we have been able to show that the soft breaking of the BRST symmetry exhibited by the 
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Gribov-Zwanziger theory can be converted into a linear breaking pj, namely the Gribov-Zwanziger ac- 
tion admits an equivalent formulation for which the BRST symmetry is only linearly broken. This is 
a useful observation, due to the fact that a linear breaking is compatible with the Quantum Action 
Principl^, while preserving the nilpotency of the BRST operator. As such, the linearly broken BRST 
symmetry can be directly translated into a system of Slavnov- Taylor identities. Furthermore, similarly 
to the case of the dimension two BRST soft breaking [9], the linearly broken formulation allows us tow 
write down additional functional identities, so that a large set of Ward identities can be established. 

Here, we prove that these Ward identities ensure the multiplicative renormalizability of both the Gribov- 
Zwanziger and Refined Gribov-Zwanziger [20^ [TT] actions within the linearly broken BRST formulation. 
As we shall see, only two renormalization factors are needed, which will be identified with the renormal- 
ization factor of the gauge coupling constant g, Zg, and of the gauge field Za- We point out that the 
linearly broken formulation enables us to establish a particular Ward identity which governs the depen- 
dence of the allowed counterterms from the Gribov parameter 7^. This Ward identity provides a simple 
understanding of the nonrenormalization properties of the Gribov parameter, expressing the fact that the 
renormalization factor Z^2 is not an independent quantity, according to the relation Z^2 = Zg ^^"^Z^^^^. 

The work is organized as follows. In Sect. 2 a brief review of the Gribov-Zwanziger action within the 
linear BRST broken formulation is provided. In Sect. 3. the set of Ward identities is established. Sect. 
4 deals with the algebraic characterization of the allowed invariant counterterm. We show that it can 
be reabsorbed into the starting action by a multiplicative redefinition of the parameters and fields, thus 
establishing the renormalizability of the Gribov-Zwanziger action. In Sect. 5 the renormalizability of the 
Refined Gribov-Zwanziger action is addressed. In Sect. 6 we outline a few observations on the break- 
ing of the BRST symmetry in a confining theory, within the context of the Refined Gribov-Zwanziger 
framework. 



2 Linearly broken BRST formulation of the Gribov-Zwanziger theory 

In this section we give a short account of the BRST linearly broken formulation of the Gribov-Zwanziger 
action. Let us start with the conventional formulation of the theory, which exhibits a soft BRST breaking. 



2.1 The Gribov-Zwanziger action and its soft BRST breaking 

As already mentioned, the Gribov-Zwanziger action implements in a local way the restriction of the 
domain of integration in the functional integral to the Gribov region fi, which is the set of all transverse 
gauge configurations for which the Faddeev-Popov operator, —d^ {8^6"''' + gf°''^^A'^^), is strictly positive 

n = {A;, d,A; = 0, -d, [d,6-' + 5/"^"^^) > } . (2) 

^For an introduction to the Quantum Action Principle, see [19] and references therein. 
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The boundary, dQ, of the region 0,, where the first vanishing eigenvalue of the Faddeev-Popov operator 
shows up, is called the Gribov horizon. The action of the theory is given by the following local expression 

Sgz = j d^x{^^F-^^F-^ + ih'^d^Al + -c''d,Dfc^ 

= \ j d^xF^^F-^, + sjd^x {cfd.Al - ulfd^Df^^^) + S, , (3) 

with 

S, = j d'x (7^ g r^^Ali^'^^ - - d{N' - 1)7^) , (4) 

where is the number of colors and d = 4 is the spacetime dimension. The fields {(p1^, 93^^) are a pair of 
complex conjugate commuting fields, while {uj^^uj^^) are anticommuting fields. Each of these field has 
4(A^^ — 1)^ independent components. The Gribov parameter 7^ is not an independent parameter of the 
theory. It is determined in a self-consistent way by the gap equation |3l HI [5] 

d£, 



■-vac 



0, (5) 



where E^ac is the vacuum energy 

e-^""'^ = j [#] e-'^Gz ^ (6) 

and \d(j)\ stands for integration over all fields entering the expression ([3|). In the absence of the term S^, 
which would imply the removal of the Gribov horizon, the action ([3]) turns out to be equivalent to the 
ordinary Faddeev-Popov action [5]. As such, it displays exact BRST symmetry 

sAl = -D'^^c' = -{8^6'^' + gr'A'^y'' , 
2 

sc" = ib" , sft'^ = , 

with 

s I d'x + . (c'^d.Ai - v'ff ) ) = . (8) 

However, in the presence of the horizon, thus 7^ 7^ 0, the Gribov-Zwanziger action is not left invariant 
by the BRST transformations, eqs.dT]), which are broken by the term 5^ 

sSgz = sS, = 7' 1 d'x {-g r'^{Dfd'){^l^ - ^'p+gr'^^AlJ^) . (9) 
Notice that, being of dimension two in the fields, the breaking term is soft. 
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2.2 Converting the soft breaking into a linear breaking 

As it has been pointed out in [1], the soft breaking of the BRST symmetry exhibited by the Gribov- 
Zwanziger action can be converted into a linear breaking. To that purpose, we have introduced two sets 
of BRST quartets of auxiliary fields, (C^l^, A^t, ^^^^ (P^' ' ^^t' <' /^^t). transforming as 

/7afe _ \ab \ab _ r, ab _ nab nab _ p. 

sptu = >'%^ sAf. = 0, sf)% = pf,, sp% = Q. (10) 
Here, {Xfu,vfv^^fu^vfv) are commuting fields, while {C^l,Cf^, p^i, p^^i) are anticommuting ones. 

Thus, as discussed in details in \X\, the Gribov-Zwanziger action ([3]) can be cast into the equivalent 
form 

+ 1 d^x + i^~\fX'5,,) . (11) 

To prove the equivalence between the two formulations, ^ and (fTTj) . one proceeds [T] by integrating out 
the fields (AJ^^, r/°^) and (A^^, r/^^), a task easily done due to the way in which these fields enter the action 
(fTTj) . Further, the fields (C°^,C^[^), (p^t; /^/ilt) "^an be completely decoupled from the theory by making 
a suitable field redefinition [1]. This leads to the equivalence between the two formulations eq.(l3]), and 
eq.dll]) . 

By making use of the BRST transformations ([7]) and (jlOp . it is apparent that, instead of a soft breaking, 
the action (jlip exhibits a linear breaking of the BRST symmetry, i.e. the resulting breaking term is 
linear in the fields, namely 

sS'tz = l' Jd'xd'^'d.X^l. (12) 
Notice also that the BRST operator s, as defined by eqs.jT]) and eqs. ffTO]) . is nilpotent 

= . (13) 

As a linear breaking is compatible with the Quantum Action Principle [19], it can be directly translated 
into Slavnov- Taylor identities. This is the task of the next section. 



3 Ward Identities 

In order to derive the Ward Identities obeyed by of the action (fTTI) . it turns out to be useful to introduce 
the multi-index notation [3 H El El [9] . 

,^0. , ,a -a\ 

(€"■■ A"- r?"- C"- ) 

V^/^«' V*' 'M*' M*/ ' 



,afe -ab , .ab 



-ab 



^ab ^ab 



-ab 



ab pa 



ravi '^au) 'I, 



ab -ab ab 
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where i = {a,fi} = with / = d{N'^ — 1) and d = A being the dimension of the Euchdean 

spacetime. The expHcit expression of the action Sq^ reads 



nlin 



1 



(15) 

Following the set up of the Algebraic Renormalization |19] . in order to account for the nonlinear BRST 
transformations of the gauge field A'^ and of the ghost field c", we introduce a pair of BRST external 
sources (ilJ^,L") coupled to the composite operators {sA^^) and (sc"), namely 



5'cxt 

Thus, the complete action 

fulfills the following set of Ward identities: 



olin I Q 
^GZ ~r *-'ext 



(16) 
(17) 



The linearly broken Slavnov- Taylor identity: 



5(S) = 7^ / d'x 5'^' V'JLC , (18) 



where 



As one easily figures out, the presence of the linear breaking term in (jlSp is inherited from the 
linearly broken BRST invariance, eq. (|12p . Also, as a consequence of the nilpotency of the BRST 
operator s, (I13p . it follows that the linearized operator By, defined as 

J \6Al^6n<^ Sn<^SA<j^ dL'^dc^ Sd^ (5c« 

ByBy = , (21) 

is nilpotent too. As we shall see later, this operator will enter the algebraic characterization of the 
most general invariant counterterm. 
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The equations of motion of the fields 6", c", C^^, /O^j, 7/|Jj and A^^: 



§ + ^^ir-o, (22) 



= -A;1. + I'S'^'S.Jl , ^ = -r,i + ^'S^%u5l , (23) 



— ^ = /d^xJ'^V^L^^ + A^.). (25) 



The parametric Ward identity for the Gribov parameter 7: 

Notice that the right hand side of eq. (j25p is hnear in the fields. Again, this term is a linear 
breaking, not affected by the quantum corrections. As already mentioned, this Ward identity has 
a special meaning. It implies that the allowed invariant counterterm does not depend explicitly 
on the parameter 7^. As such, the Ward identity (j25p provides a very simple understanding of 
the nonrenormalization properties of the Gribov parameter 7^, as encoded in the relationship 
Z^2 = Zg ^^'^ Zj^^^ . We also point out that the possibility of writing down the Ward identity (f25]) 
is a byproduct of the linear breaking formulation of the Gribov- Zwanziger action. 

The local equation of motion of (p^: 



where A^* is a linear breaking, given by 



= -av^ - - d,nu - grr'A^si^ . (27) 



• The local equation of motion of w": 
where 

As^ = aV + 9^c^ + a^p«, . (29) 



The local equation of motion of 99?: 



with 



A«^ = -d'^^ + d.XI,, + d,y^, - j'gr'^Al5l^ . (31) 
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The local equation of motion of w^: 
and 

A-^ = -d^of + d,p';,, + d,c;,. (33) 

Notice that all breaking terms in eqs. ([29]) . (f3T]) . (j33]) are linear. 

The integrated Ward identities involving the auxiliary and the Faddeev-Popov ghost fields: 



M'^ c«— + ^?— +7^^^^ (35) 



The linearly broken ghost Ward identity [19] : 

^'^(S) = , (36) 

where 

y j''\5c- "db'^^^^ 'm^^ ^'du'?^^^ ^^^'dpi 

and the linear breaking A° given by 

= / d4x5/"''^(l^^A^ - L^c^ - • (38) 



• The linearly broken identity involving only the auxiliary fields: 

Mi (S) = -7' / d^x 6^%,Ji, {-p% + C^) , (39) 

where 

(40) 
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The linearly broken global symmetry U{f), f = 4(A^^ — 1): 

Q., (S) = j^S-\, J d'xiSlX^, - 5l\%) , (41) 



where 





5 
SC. 


na 


5 


5 




6 


) 



c r e r 

rt" — - ^n"--— A" — l-fi" — A" ■ - 1 (42) 

The linearly broken rigid identity: 

>V'^(S) = 7^^'^ / d'x {\% + r?^^^) , (43) 

with 



-b ^ , b ^ , -b S , b ^ , -b ^ , b ^ , nb ^ 



4 Stability and invariant counterterm 

4.1 Algebraic characterization of the invariant counterterm 

In order to characterize the most general invariant counterterm which can be freely added to all orders in 
perturbation theory we follow the set up of the Algebraic Renormalization [12] and perturb the classical 
action S by adding an integrated local polynomial in the fields and sources, E'^*, with dimension bounded 
by four, and with vanishing ghost number. We demand thus that the perturbed action, (S + eS"^*), 
where e is an expansion parameter, fulfills, to the first order in e, the same Ward identities obeyed by 
the classical action S, i.e. equations (jlSp - (j43p . This requirement gives rise to the following constraints 
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for the counterterm T, 



ct. 



^sS'=* = 0, (45) 



— ^=0' T^ = 0' 1J7T = 0, (46) 

(47) 
(48) 

.^^(S^*) = 0, HS^s(S^*) = 0, Z|.(S^*) = 0, J^(S-*) = 0, (49) 

g«(S^*) = 0, Mj(S^*) = 0, Q,j(S^*) = 0, W(S^*) = 0. (50) 

Equations P6]) imply that the counterterm is independent from the fields (t*^, ^/JJlj, A^j, C|^j, p^J and from 
the parameter 7^. From equation (j47p it follows that c° and 0^ enter only through the combination 
{d^cf + r2^). We also notice that the functional operators in equations (09]) are those associated with the 
Ward identities ([30]), ([MD and ([35]), namely 





0, 










+ 9^ 








H^'(S^*) = 


0, 


HSnS^*) = 0, 








= 0, 


Z|.(S^*) = 0. 


, j^(s-*) 


= 0, 


M,-(S^*) = 


0, 


Q.j(s^*) = o, 




= 0. 



_„fabc_^J__ fabc_^J_ 



Further, the trace of the operator Qij in eq. (f50]l . defines a new quantum number, the Q-charge, displayed 
in table [H Finally, reminding that the operator By. is nilpotent, eq. (l2ip . from eq. ()45p it follows that E'^* 
can be characterized by looking at the cohomology of the operator By. in the space of the integrated local 
polynomials with dimension four and zero ghost number and Q-charge. Thus, according to the general 
results on the cohomology of Yang-Mills theories [19], we start by writing the most general invariant 
counterterm as 



vet _ oo 



j d^xF-,F-, + BYA(-^K (55) 



The first term in the r.h.s of equation (j55p represents the nontrivial part of the cohomology of By, with 
ao being an arbitrary coefficient, while A^"^) is an integrated polynomial in the fields and sources with 
dimension 4, ghost number —1 and vanishing Q-charge. 
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A 


h 


c 


c 


VL 


L 


95 




UJ 




A 


A 


V 


7) 


C 


C 




P 


dimension 


1 


2 





2 


3 


4 


1 


1 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


ghost number 








1 


-1 


-1 


-2 








1 


-1 














1 


-1 


1 


-1 


Q-charge 




















1 


-1 


1 


-1 


-1 


-1 


1 


1 


1 


-1 


1 


-1 



Table 1: The quantum numbers of fields and sources of the theory 



By using the equations (|i6]) and (jl7|) , we find 

^'^ + Pt'"' + f^t' u^tut^u^ + iotu^'Muf) , (56) 

where {k = 1,...,10) are constant arbitrary coefficients, while Pf'"^'^ (/ = 1,2,3,4) stands for an 
invariant tensor of rank 4. Using the remaining constrains on E*^*, and taking into account the following 
useful relations: 

=-SS*s, {^s,HS^}=H-^, [Bj:,Ih]=-Jh, {Ss,^"} = W^ (57) 
it turns out that A^"^^ depends only from one free parameter ai, namely 

A(-^) = ai 1 d'x ( {d,c^ + n;) A^ + Cjtd,Df^\ + n%Dfu^\ - C%D'^^^\ + C^,f?^, ) . (58) 

Therefore, for the final form of the most general allowed invariant counterterm one finds 

-^tdi^ti + ^'idAi + '^"^V? - . (59) 

4.2 Renormalization and Z-factors 

It remains to check that the counterterm (j59p can be reabsorbed in the starting action S by means of a 
multiplicative redefinition of the fields, sources and parameters of the theory, i.e. 

2(0, $) + e S^*(</., $) = S(,^o, ^-o) + O(e') (60) 

where (0, 0o) denote the renormalized and bare fields, respectively, while ($, $0) stand for the renormal- 
ized and bare sources and parameters, i.e. <^ = (L", fi^, (7, 7^). The renormalized and bare quantities are 
related each other as 

</)o=^y'0, ^Q = z^^, (61) 
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where the Z's are renormahzation factors. One can easily proove that eq. ()60p is in fact fulfihed by writing 

1/2 

the bare action in terms of two independent renormahzation factors Zg and Z , i.e.: 

S(<Ao,^o) = 1 d^x{^{d^,Al)d^Al^ZgZT9r''\d^,Al)A\Al^\zlZ^^^^ 

+iz'J%%,A';, + z-^z^^/ Va^c'^ + gr''{d^c^)A<;,c' 

+Zg'z~'/'x^^,d,^t + 9^''>^'',^A';,^^ - z-^zi"^niid,^i 

-grXiA^^'i + Zg^Z~^''gr'%,{d,c'^)oj\ + g^r"^nPf)l,A^d>oj\ 



-Zg^Z'^"\l^l^ - Zr^^Z-^I^Cl^^^l - gr'Cl.A^^Ji - Z-^I^Z-^^pl^^^l 

-gr\%A'^^\ + + z-^zi^i^ci^p^^, + Z-^^Z-2l\r'%Ad^c'^)H}'i 

9 
2' 



XrK^ - Zg'Z^'/'n-d.c'^ - gn-r'A'^^c^ + ^r'^l'^c'A . (62) 



The two independent renormahzation factors Zg and Z^'^ are related to the coefficients ag and oi in the 
following way: 

Zq = 1 — -ao , 



ZT = l + .(? + ai). (63) 

(64) 



2 

2 

1/2 

All remaining renormahzation factors can be expressed in terms of Z^ and Zg as 

Zb = Z ^ , = = Z^ = Z^ = Z\ = Zf^ = Zg ^Z ^ , 

^1/2 _ ^1/2 _ ry-\ ryl/2 _ ^1/2 _ ^"1/2 ^ _ ^7-1/2 ^-1/4 ^^1/2 _ ^7 

^Cb ~ ~ 9 ' P ~ ~ ^ A ' ~ ^9 ^ A ' ~ ^9 ' 

■2^L = -^p^ = Z , = = ZgZ , = Zg ^ Z ^ . (65) 

In particular, we point out the two relations: Z^Z^J^ Zg = 1 and Z^2 = Zg ^^'^Zj^^'^. The first one 
expresses the nonrenormalization theorem of the ghost-antighost-gluon vertex. The second relation tells 
us that the renormahzation factor Z^2 is not an independent parameter of the theory. Within the linear 
breaking formulation, this result is a consequence of the Ward identity ()25p . This ends the proof of the 
multiplicative renormalizability of the Gribov-Zwanziger action in the linear breaking formulation 



5 Generalization to the Refined Gribov-Zwanziger 

The results established in the previous section can be generalized to the so-called Refined Gribov- 
Zwanziger action [201 [H], which takes into account the effects of the dimension two condensates (^^^^) 
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and (ifl^if'^ — uj^uj^). The presence of these condensates can be encoded in the starting action by modifying 
the Gribov-Zwanziger action as 

Srgz = Sgz + j d^x {^f^f - Cofiof) + ^ J d^x A^A^ , (66) 

where the mass parameters /i, m are dynamical parameters related to the aforementioned dimension 
two condensates, see [22] for an updated discussion. The Refined Gribov-Zwanziger action gives rise 
to a gluon propagator which remains suppressed in IR region. However, unlike the propagator of the 
Gribov-Zwanziger theory, eq.([l]), it attains a non- vanishing value at the origin in momentum space 

Al{k)Ati-k)) _ = , , l!/ f V - ^) , (67) 



RGZ k"^ + (m^ + ^2)/j2 _j_ ^2^2 _^ ry4: y 1^" ^2 

Expression (j67p in in agreement with the recent lattice simulations [231 [2^ [25l [26] . A similar behavior 
for the gluon propagator has emerged from the studies of the Schwinger-Dyson equations [271 HSj , see 
also p9j for the introduction of the concept of the dynamical gluon mass. 



In order to discuss the Refined Gribov-Zwanziger action, eq. (|66p . within the linearly broken BRST formu- 
lation, we analyze the two terms, {(Pi^l — ^^fwf) and ^^^4^, separately. Let us start with (<^fv3° — tDfwf), 
which is BRST invariant. In fact 

d^x - Ofuf) = sjd^x i^^uf) . (68) 

This term will thus not affect the BRST breaking. Moreover, it is easy to check that all Ward identities 
in eqs.([l8]), ([22]), ([23]), ([25]), ([26]), ([28]), ([30]), §21, §6^, ([39]), (BII), & remain valid, up to possible 
irrelevant linear breaking terms. Consider, for example, the equation of motion of (^°, eq. ()26p . With 
the inclusion of {'f'^^i — Cofu'^), the breaking term A^*, eq. fpT]) . gets modified by the addition of a linear 
term, i.e. 

A^*^A| + ;uV- (69) 

As such, the Ward identity (j26p retains its full validity. A similar argument applies to the other identities. 
Of the whole set of Ward identities established in the previous sections, only the two Ward identities in 
eqs. (j3l]) and ([35]) cannot be maintained due to the introduction of the term {(pffi — Ofujf). Nevertheless, 
it turns out that the lack of these two Ward identities does not prevent us to prove the renormalizability 
of the model. 

The second term, A^A^, in expression (j66p is more delicate, since it is not BRST invariant 

s J d'x \aIAI = - fd'x Ald.c'^ . (70) 

Moreover, as the induced breaking is a soft breaking, we can repeat the previous argument and convert 
this breaking into a linear one. Let us briefly sketch how this is done. We introduce a BRST quartet of 
auxiliary fields (£,£,a, (3), where are a pair of anticommuting fields, while {a, (3) are commuting, 

s£ = a , 

sa = , 

= £, 

s£ = 0. (71) 
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Thus, the term / d'^x ^JJ^J^ is replaced by the equivalent term Sm 
Sm = s J d^x Af^A"^ - +n? j 



/ 



d^x [ \aAlAl + SA^^df.c" + ££ - /3{a - m') ] . (72) 



Notice in fact that from the equation of motion of /3 we have: 

5Sm 



513 

Therefore, making the substitution a = rm?, it follows 



(a - m^) = . (73) 



2 

m 



Sm = j d'x [j-AlA^ + £{A1 d.c'^ + £)J . (74) 

Now, performing a linear shift in the £ variable, with unity Jacobian, 

£^£-A1d^c\ (75) 

we obtain ^ 

S^^ j d^x (^A^Al^ + j d^x££. (76) 

As the last term, / d'^x££, is completely decoupled from the rest, the fields {£,£) are harmless and can 
be integrated out, so that the starting term tit? J d'^x A'^^A'^ is recovered. The soft breaking in eq.(|70p 
has thus been converted into a linear breaking, as it is apparent from 

sSm = rn^J d^x£{x). (77) 

The linear breaking ([77|) can be now easily incorporated into the Ward identities ([H]), ([22]) . ([23]) . (p^ . 
(I25I1 . (1261) . aiidihe previous algebraic analysis can be repeated. We 

limit here only to state the final result, confirming that the refined action (j66p is renormalizable within the 
linear breaking formulation. In particular, the nonrenormalization theorem of the ghost-antighost-gluon 
vertex, ZgZcZ^j^ = 1, and the relation = Zg ^^'^Zj^^^ remain valid. 



6 Considerations on the breaking of the BRST symmetry, Gribov 
horizon and dimension two condensates 

In this section we provide a few remarks on the breaking of the BRST symmetry within the context 
of the Refined Gribov- Zwanziger framework, obtained by taking into account the nonperturbative ef- 
fects of the Gribov horizon and of the dimension two condensates (^d^^dj^), {(p^'^'l — wfu;"). As we shall 
see, the Refined Gribov-Zwanziger set up enables us to argue that the breaking of the BRST symmetry 
might not be in conflict with the confining character of the theory. Rather, it allows us to put forward 
the idea that the nonperturbative effects of the Gribov horizon and of the dimension two condensates 
conspire in such a way that the resulting BRST breaking does not spoil at all the consistency of the theory. 
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Let us start by reminding that, in the Refined Gribov-Zwanziger framework, the gluon propagator, 
eq. ()67p . violates positivity while attaining a nonvanishing value at the origin [20^ lllj. Concerning the 
ghost propagator, it essentially keeps a free behavior at very low energies [201 [H], namely 



As already mentioned, this behavior of the gluon and ghost propagators is in good agreement with the 
most recent lattice data [231 IMl 123 |26]. In particular, expression ([67|) provides an accurate fit of the 
lattice data up to 1.5 GeV [26]. To some extent, expression ([78]) can be taken as evidence of the fact 
ghosts behave freely at very low energies, i.e. they decouple in the deep infrared. 

This feature is corroborated by the observation that, in the Refined Gribov-Zwanziger framework, the 
renormalization group invariant ghost-antighost-gluon effective coupling, cngh-gi 

where and J are the gluon and ghost form factor^ [30j, goes to zero in the deep infrared, i.e. 
ctgh-gi{k'^)\^^Q 0. The coupling agh-gi{k'^) is directly related to the ghost-antighost-gluon vertex of 
the Gribov-Zwanziger action and of its refined version. The fact that it goes to zero at very low energies 
can be seen as a further signal that ghosts get less and less interacting in the deep infrared. It is worth to 
emphasize that we are not claiming here that the quantity agh^gi{k'^), as defined above, has the meaning 
of the physical coupling constant of Yang-Mills theory. To our knowledge, the problem of the physical 
definition of the YM coupling constant is far from being solved and is currently under intensive debate. 
See [31] for a recent proposal. 

We are now ready to present our considerations about the breaking of the BRST symmetry within 
the Refined Gribov-Zwanziger framework. 

6.1 Renormalizability 

The first issue which we point out is that of the renormalizability. One might argue that the breaking of 
the BRST invariance could jeopardize the renormalizability of the theory. Though, as we have discussed 
in the previous sections, this is not the case. The breaking of the BRST symmetry induced by the Gribov 
horizon and by the dimension two condensates is fully compatible with the renormalizability of the theory. 
This is well captured by the present linear formulation in which the linearly broken BRST invariance can 
be directly translated into useful Slavnov- Taylor identities, eq. ijlSp . while keeping the nilpotency of the 
BRST operator s, eq. (fT3]l . 

6.2 Gluon confinement 

In a nonconfining gauge theory the characterization of the physical subspace heavily relies on the notion 
of the BRST charge, Qbrst, and thus on the exact BRST invariance of the theory. We remind here that 

■^We remind here that the possibihty of introducing the RGE invariant quantity oigh-giik^) ~ V{k'^)J'^{k^) rehes on the 
nonperturbative vahdity of the nonrenormahzation theorem of the ghost-antighost-gluon vertex, ZcZgZ]/'^ = 1, see ref. [30] 
for a lattice investigation. 

^The gluon and ghost form factors, Vik^), J''^{k^), are introduced through the expressions 



c''ik)c\-k))nGZ 



(78) 





■ab 



(79) 
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the BRST charge Qbrst enables us to remove from the asymptotic physical subspace the unphysical 
polarizations of the gauge field, i.e. the longitudinal and the temporal modes, as well as the negative 
norm states associated to ghost and antighost states. The resulting subspace, identified by means of the 
cohomology of the BRST charge Qbrst in Fock space, contains only positive norm states corresponding 
to the physical transverse polarization of the gauge field. All this is well understood and applies to non- 
confining theories for which the asymptotic fields can be safely introduced, having a direct relationship 
with the excitations of the physical spectrum. 

In a confining theory things are completely different. Gluons are not part of the physical spectrum, 
due to color confinement. One needs thus a mechanism which enables us to remove them from the phys- 
ical spectrum. Here, we argue that such mechanism is provided by the combined effects of the Gribov 
horizon and of the dimension two condensates. These effects conspire in such a way that gluons become 
unphysical in the low energy region, as it follows from the expression of the gluon propagator (j67p which 
exhibit positivity violation, while displaying complex poles [20\ lllj. As a consequence, transverse gluons 
cannot correspond to physical excitations of the spectrum. On the other hand, in the Refined Gribov- 
Zwanziger framework, ghosts seem to decouple at very low energies. It should be pointed out that the free 
behavior of the ghosts in the deep infrared is essentially due to the dimension two condensates. Neglecting 
these condensates, would lead to a very different behavior of the ghosts in the infrared. In fact, the ghost 
propagator would be enhanced in the deep infrared, {c°'{k)c^{—k))Gz\k^o ~ while the ghost-antighost- 
gluon effective coupling Ugh-giik"^) would attain a nonvanishing value, i.e. agh-gi{k'^)\i^^Q ad^^ [32j. 
All this indicates that, without the inclusion of the dimension two condensates, ghosts would be strongly 
interacting in the infrared. In summary, in the Refined Gribov-Zwanziger framework, the Gribov horizon 
and the dimension two condensates provide us a mechanism for the removal of the transverse gluons 
from the spectrum, while turning ghosts almost free at very low energies. In this sense, we would have 
a tendency to believe that the breaking of the BRST symmetry does not compromise the consistency of 
the theory. 

6.3 The physical spectrum 

We come now to the pivotal issue to be addressed by any framework aiming at facing a confining theory: 
the analytic characterization of the physical spectrum. Even if a possible mechanism for the removal 
of the gluons from the spectrum of the theory has been worked out, this has to be regarded only as a 
preliminary step. The real challenge is the physical spectrum of the theory and the analytic calculation 
of the masses of the physical excitations which, in the present quarkless case, are glueballs |33] . 

Prom a theoretical and first principle point of view, the physical spectrum of the theory can be accessed 
through the evaluation of the correlation functions, {0{k)0{—k)), of a suitable set of local composite op- 
erators {O}. This is a highly nontrivial task, given that the correlator {0{k)0{—k)) has to be evaluated 
with a confining, positivity violating, gluon propagator which exhibit complex poles, see expression (j67p . 

On physical grounds, we expect that, in order to be a good candidate, a local operator O should dis- 
play several features. Of course, we cannot be exhaustive here. We quote, for example, the following 
properties: 

• i) the operator O should be a colorless operator and we should be able to consistently evaluate the 
correlation function {0{k)0{—k)) . We look thus at renormalizable operators. 
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• a) the correlation function should satisfy the Kallen-Lehmann spectral representation 



poo 

{0{k)0{-k)) = dT 
Jo 




(80) 



with a positive spectral density /o(r). 

Condition i) is a basic requirement. Certainly, the BRST plays a very important role here. We remind 
in fact that, even in the presence of the breaking, the operator s keeps its nilpotency, = 0, eq. (fT3]) . so 
that one still has the notion of the BRST cohomology. Moreover, looking at the BRST transformations, 



eqs.dTD, m, one reahzes that all fields (c^ 6", a;^^ ^f, iff) and {Cf,, \%X'u. C% 1^, nfu^pfu) 



give rise to BRST doublets [19]. As a consequence, the cohomology classes of the BRST operator with 
zero ghost number are given by gauge invariant colorless operators built with the field strength F^^^ and 
its covariant derivative [19] • Although we lack exact BRST invariance, there is an important observation 
which can be made about the local composite operators belonging to the cohomology of s. Remarkably, 
these operators can be proven to be renormalizable, and this even in the presence of the BRST breaking, 
see [18] for a detailed discussion of the renormalizability of the operator F'^ix) = F^^{x)F^^{x). In other 
words, the existence of a BRST breaking does not invalidate the renormalizability of the gauge invariant 
operators belonging to the cohomology of the BRST operator. This feature can be taken as a support in 
favor of the fact that operators belonging to the BRST cohomology might be useful for the investigation 
of the spectrum of the theory [SU [35l [36]. Evidently, due to the presence of the breaking, the BRST 
exact piece of the renormalized version Oren of a given operator O has to be treated with due care when 
computing the renormalized correlator {Oren{k)Oren{—k)) ■ 

Nevertheless, to our understanding, the requirement of renormalizability represents only a first step 
towards the identification of a meaningful composite operator. The fulfillment of the spectral represen- 
tation, eq. llSOp . with a positive spectral density seems to be a key ingredient in order to provide a link 
with the physical spectrum. In fact, the knowledge of the spectral density 



would give us access to the mass M of the excitation associated to the composite operator O and to the 
threshold tq for the production of multiparticles states.. In addition, the spectral representation ([80]) 
will ensure that the correlation function {0{k)0{—k)) can be Wick rotated from the Euclidean to the 
Minkowski space, a feature of fundamental importance for a correct physical interpretation of the theory. 
Though, till now, the existence of the spectral representation seems to be uncorrelated with the issue of 
the BRST symmetry. In particular, it is not easy to figure out a possible relationship between BRST and 
the positivity of the spectral density. We remind here the fact that the correlation function {0{k)0{—k)) 
has to be evaluated with a confining gluon propagator which violates positivity. From this, one can really 
appreciate the hard task of accessing the spectrum of the glueballs. In this context, it is worth emphasizing 
that the analytic structure emerging from the use of the refined gluon propagator (I67p looks quite useful. 
As discussed in [10], the complex conjugate poles associated to expression (i67|) combine in such a way 
that the final result does display a nice spectral representation with positive spectral densities, at least 
for a certain class of composite operators. This feature has enabled us to extract numerical estimates for 
the masses of the three lightest glueballs, 0^^,2^^,0 ^ [36], which turn out to be in good agreement 
with the available lattice data, see [33]. This is an encouraging result, confirming that the interplay 
between the Gribov horizon and the dimension two condensates might be helpful in order to access the 
glueball spectrum. Though, we have the impression that, in a confining theory, BRST alone would not 



p(t) 



7^ 6{t - M^) + e{T - ro)p(r) 




(81) 
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be sufficient for a complete characterization of the spectrum. Some additional requirements gauranteeing 
the validity of the spectral representation and of the positivity of the corresponding spectral density seem 
to be needed, especially in view of the fact that explicit calculations have to be carried out with a gluon 
propagator which violates positivity. 



7 Conclusion 

In this work, the soft breaking of the BRST symmetry due to the Gribov horizon and to the dimension 
two condensates has been investigated. This breaking can be converted into a linear breaking, a feature 
which has enabled us to write down a useful set of Ward identities. These Ward identities ensure the 
multiplicative renormalizability of both Gribov- Zwanziger action and its refined version. 

We underline that the behavior of the gluon and ghost propagators obtained from the Refined Gribov- 
Zwanziger framework turns out to be in good agreement with the most recent lattice data. Recently, 
the refined gluon propagator has been employed to estimate the masses of the three lightest glueballs, 
0++,2++,0"+. Also here, the results are in agreement with the available data. 

Although much work is certainly needed in order to achieve a satisfactory understanding of the role 
played by the BRST symmetry in a confining theory, the encouraging results obtained so far within the 
Refined Gribov-Zwanziger framework have enabled us to put forward a few observations in favor of the 
idea that the breaking of the BRST symmetry, in both soft and linear formulations, does not jeopardize 
the consistency of the theory. 
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